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ABSTRACT 

exact  solution  fur  ray  paths  nestr  a  sound  channel  with  parabolic 
velocity  profile  Is  given  for  rays  crossing  the  sound  channel  axis  at 
small  angles.  The  convergent  zone  for  such  rays  for  a  source  on  axis 
Is  compared  with  the  numerically  computed  results  of  M.  A.  Pederson, 
J.  Acoust.  Soc>  Am.  ,  ^3,  465,  (1961).  Two  simple  straight  line  velocr 
Ity  gradient  approximations  are  compared  with  the  parabolic  case  for 
the  purpose  of  observing  their  accuracy.  Both  cases  show  surprisingly 
good  agreement  over  a  range  of  choices  for  the  gradients. 
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I.  INTRODUCTION 


The  problem  of  computing  ray  paths  for  sound  rays  \mder  actual  oceano¬ 
graphic  conditions  is  known  to  be  difficult  due  to  the  complicated  depth  profiles 
generally  fo\md.  This  is  especially  true  imder  general  summer  conditions  in 

which  there  exists  a  strong  velocity  gradient  just  below  the  surface  channel^ . 
The  usual  approach  to  ray  tracing  is  one  in  which  the  velocity  profile  is 
approximated  by  straight  line  segments  used  over  various  depths.  Such  an 
approximation  will  obviously  lead  to  errors  in  both  the  ray  paths  and  in  the 

2 

predicted  convergent  zones  .  It  is  the  purpose  of  this  report  to  present 
a  simple  solution  given  in  Section  II,  for  the  ray  paths  and  convergent  zones 
for  sound  velocity  profiles  exhibiting  minima-  The  solution  given  is  valid  for 
rays  leaving  the  source  at  small  angles  from  the  channel  axis*  For  conven¬ 
ience  this  paper  is  limited  to  soimd  sources  on  the  sound  channel  axis-  The 
extension  of  the  source  to  off  axis  position  is  straightforward  as  long  as  the 
small  angle  requirement  mentioned  above  is  met*  Since  the  constant  gradient 
approach  to  velocity  profiles  is  so  useful  in  many  problems  in  Section  III  two 
very  simple  constant  gradient  examples  are  fitted  to  the  parabolic  profile  and 
their  convergent  zones  are  compared  with  the  parabolic  results.  The  possible 
use  of  the  parabolic  result  is  discussed  in  Section  IV,  as  well  as  the  merit  of 
the  straight  line  method.  The  possibility  of  considering  the  types  of  profiles 
by  the  above  method  is  mentioned. 


C.  B.  Officer.  Introduction  to  the  Theory  of  Sound  Transmission, 
p.  150.  McGraw-Hill  Book  Co.  New  York-  (1958). 

Melvin  A.  Pederson.  J.  Acoust.  Soc.  Am.  33,  465  (1961). 
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II-  THE  RAY  PATH  FOR  A  PARABOLIC  VELOCITY  PROFILE 

There  exists  a  very  simple  solution  for  the  ray  equation  for  small  angles 
If  the  velocity  has  a  depth  dependence  given  by 

C(Z)  =  C  (1  +ot^Z^)  .  (1) 

o 

Where  the  axis  Z  =  O  (corresponding  to  what  shall  be  called  the  X  axis)  Is 

taken  at  the  minimum.  Typical  numbers  given  by  Pederson^  Indicate  that 
for  veilues  of  reasonably  close  to  the  axis  (Z  =  O)  the  error  in  approximating 
the  velocity  profile  by  this  form  Is.  less  than  one  per  cent-  This  Is  seen  from 

the  equation  given  by  him  on  page  471  of  his  article.  From  Snell's  Law 


_ 1 _  „  cos  ftp  (2) 

Vl+  Z'^'  C(Z)  ° 

where  S  is  the  angle  of  emergence  of  a  ray  from  the  source  and  Z'“  -r^  , 

2 

It  may  be  seen  that  If  C(Z)  only  varies  slightly  from  C^  then  Z'  is  also  small. 

The  condition  on  C(Z)  is  certainly  the  oceanographic  case.  If  0^  is  small,  the 

substitution  of  Eq-  (1),  expansion  of  Eq.  (2),  suid  omission  of  higher  order  terms 
leads  to  the  differential  equation 


+  oc2z2 


(3) 


M.  A.  Pederson-  J.  Acoust.  Soc.  Am.  33,  465  (1961). 
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This  has  the  same  form  as  the  Hamiltonian  for  the  harmonic  oscillator  in 
classical  mechanics,  the  solution  of  which  is 

0 

Z  =  — ^  sin  Vi  o<r  -  (4) 

Vi 

It  is  then  immediately  obvious  that  all  rays,  independent  of  8^,  cross  the 
axis  for  the  first  time  at  the  range 

R  =  (5) 

^2  cK 

2 

This  is  correct  up  to  the  order  6^  .  From  the  numbers  given  in  Pederson's 

article  on  page  471^,  the  predicted  convergent  zone  using  Eq.  (5)  is 
R  «  5,623  yards,  which  is  in  very  good  agreement  with  Pederson's  value, 

5,539  yards.  The  discrepancy  is  due  mostly  to  the  approximation  of 
Pederson's  velocity  equation  by  a  parabolic  fimction. 


M.  A.  Pederson. 


J.  Acoust.  Soc.  Am.  33,  465  (1961). 


NATO  UNCLASSIFIED 


NATO  UNCLASSIFIED 


in.  APPROXIMATIONS  TO  THE  ABOVE  RESULTS 
USING  STRAIGHT  LINE  SEGMENTS 


With  the  above  result  giving  a  simple  range  formula  it  is  of  Interest  to 
compare  with  it  the  results  of  the  straight  line  segment  method  of  approx¬ 
imating  velocity  profiles.  The  two  different  cases  considered  are  shown 
In  Fig.  1  ■  The  first  is  a  line  segment  of  constaxit  velocity  joined  with  a  line 
segment  tangent  to  a  point  on  the  parabola.  Very  simple  calculations  given 
In  Appendix  A  show  for  this  case  a  convergent  zone 

R,  -  ^  (6) 

which  Is  90.  5%  In  range  of  that  given  Eq.  (5).  This  Is  surprisingly  good 
agreement  for  such  a  simple  approximation.  The  second  approximation 
considered  is  that  given  by  two  line  segments  as  shown  In  Fig.  1.  The  point 
nearest  the  axis  was  taken  as  the  fraction  S  of  the  height  of  the  second 
point.  The  calculated  convergent  zone  as  a  function  of  as  derived  in 
Appendix  A  is 

R,  .  ^ 

•‘(I  +  5  ) 

This  function  R2(  $  )  monotonically  Increases  over  o  A  ^  ^  1  .  It  rapidly 

rises  to  about  90%  of  the  value  given  by  the  parabolic  solution  (Eq.(  5))  and 
increases  to  about  110%  of  that  value  at  ^  -1.0.  At  the  value  3  ■  0.  61. 

R^  is  equal  to  the  parabolic  value  K/V?  04.  Thus  as  long  as  3  greater 

than  0.  4  the  agreement  with  the  parabolic  convergent  zone  is  very  good. 
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Fig.  1 

Velocity  Profiles  used  to  Approximate  Parabolic 

2  2 

Depth  Dependent  Curve  C  «  C  (1  -f  OC  Z  ). 
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IV.  DISCUSSION 


The  small  angle  approximation  has  allowed  a  very  simple  solution  to 
the  ray  path.  Thus  in  regions  where  there  exists  a  minimum  in  velocity 
profiles,  there  exists  a  region  over  which  this  solution  may  be  used.  The 
limitations  imposed  by  actual  cases  would  be  due  to  the  actual  width  of  the 
"parabolic"  region  in  the  sound  channel.  If  the  deviation  from  parabolic 
shape  is  too  great  over  a  small  width  then  in  order  to  include  a  sufficient 
range  of  ray  angles  the  method  would  have  to  be  extended  by  connecting 
other  types  of  gradients  to  the  peirabolic  one.  It  should  be  noted  that  when 
considering  ray  paths  due  to  the  sound  channel,  it  is  only  necessary  to 
consider  the  velocity  profile  on  one  side  of  the  axis  since  the  two  sides 
may  be  treated  separately.  Thus  much  of  the  actual  asymmetry  of  the 
profile  may  be  Ignored. 

The  approximation  methods  used  in  Section  lU  give  remarkably  good 
agreement  with  the  parabolic  results  of  Section  II.  This  agreement  makes 
one  more  confident  that  simple  models  of  complicated  profiles  are  still 
useful  for  some  purposes.  Pederson's  comment  concerning  anomalous 
caustics  is  not  meaningful  for  the  two  velocity  gradient  case.  It  does  be¬ 
come  important,  however,  if  one  or  more  gradients  are  added  to  the  cases 
considered  here.  One  additional  gradient  added  to  case  one  in  Section  III 
produces  another  convergent  zone  which  was  fotmd  to  wander  about 
drastically,  depending  on  the  relative  depths  of  the  layers.  The  same 
convergent  zone  given  by  Eq.  (6)  is  still  present.  Thus,  as  Pederson 
pointed  out,  caution  must  be  used  in  interpreting  the  results. 

I 

It  is  interesting  to  note  that  for  a  velocity  versus  depth  function  of  the 
form 


C(Z)  -  C  (1  +  f  (Z)  ) 
o 


(8) 
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where  f(Z)«l  for  all  Z  the  differential  equation 

i  +  «z)  -  i  «2„) 

for  rays  is  obtained  from  Snell's  Law  for  small  angles.  Z^  is  the  depth  of 

the  source  and  0^  is  the  angle  with  the  horizontal  of  the  outgoing  ray.  This 

equation  has  the  form  of  a  Hamiltonian  for  arbitrary  potentials  Thus 

many  of  the  presently  used  solutions  of  classified  mechanics  may  be  useful 
in  ray  tracing  problems  (for  small  angles).  One  of  the  authors  (R.  R.  Goodman) 
is  currently  Investigating  this  possibility. 
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APPENDIX  A 


In  this  section  the  convergent  zones  will  be  calculated  for  the  two 

gradient  cases  shown  in  Fig.  1-  The  source  is  considered  on  the  sound 

channel  axis.  First  the  distance  at  which  a  ray  leaving  the  source  at 

angle  6  crosses  the  axis  must  be  found-  Call  this  distance  R(9  ). 
o  j  o 

The  result  for  case  one  is  easily  found  by  well  known  means  to  be 


"■'“o’ 


2  d 
tan  0 

o 


2  C 

o 


tan  0 

o 


(A-1) 


where  d  is  the  height  of  the  constant  velocity  region  and  g^  is  the 
velocity  gradient  above.  Since  g^  is  the  tangent  to  the  slope  of  the 
parabola  at  some  point  Z  =  h-  By  Eq.  (1)  it  is  seen  that 

gj  =  .  (A-2) 

The  value  d  is  also  easily  seen  to  be  h/  2.  In  order  to  calculate  the 
position  of  the  convergent  zone  it  is  only  necessary  to  find  the  value  of 

dRi 

0  such  that  -tti”  =  O  •  This  value  is  seen  to  be 
o  d  0 

o 


tan  0 

o 


(A-3) 


C.  B.  Officer.  Introduction  to  the  Theory  of  Sound  Transmission , 
Chapter  II.  McGraw-Hill  Book  Co.  New  York.  (1958). 
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When  Eqs-  (A-2),  (A-3)  and  d  =  —  are  substituted  into  (A-1)  the  convergent 
range  is  seen  to  be 


R 


1 


(A-4) 


The  second  case  is  calculated  along  similar  lines-  The  range  at  which 
a  ray  crosses  the  axis  is  given  by 


2C  2C  2C 

)  =  — ”  tan  0  -  - tan  ©  + -  tan  ©  (A-5) 

2  o  gi  o  gj  g^ 

where  g  and  g  are  the  two  gradients  is  seen  in  Fig-  1-  C.  and  ©  are  the 

X  A  X 

values  of  the  velocity  and  angle  of  the  ray  respectively  at  the  depth  of  inter- 

dR2 

section  of  the  two  line  segments.  Setting  —  =  O  .  using  the  fact  that 

o 


d©  C,  sLi  © 

_  _  _ ^  _ o 

dO  C  sin  © 

o  o 


(A-6) 


which  follows  from  Snell's  Law, an  equation  for  ©^  is  obtained  which  is 


(A- 7) 
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Since  the  points  h  and  Jh  on  the  •’•elocity  profile  determine  .  g^  and 
Cj,  it  is  easily  shown  that 


gl  =  Jh 

gg  =  ^  (1  +  1  )h  (A-8) 

C  =  C  (1  +  «  ^  5  . 

1  O 

Substituting  Eqs  (A-7),  (A-8)  and  Snell's  Law  into  (A- 5),  the  range  of  the 
convergent  zone  is  found  to  be 


Tj  2  V2t  (T  +2)' 

“2  “  o<(l  +  J  ) 


(A-9) 
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